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An analytic theory is developed for the density of states oscillations in quantum wells in a magnetic 
field which is tilted with respect to the barrier planes. The main oscillations are found to come from 
the simplest one or two-bounce periodic orbits. We calculate their period and stability analytically 
and find an infinite sequence of destabilizations followed by restabilizations as the chaos parameter 
increases. This phenomenon explains the re-entrant frequency-doubling of the density of states 
peaks observed in recent magnetotunneling experiments. 

PACS numbers: 05.45. +b, 72.15.Gd, 73.20.Dx 
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Most of our intuition about the level structure of quan- 
tum systems comes from the consideration of hamiltoni- 
ans with high symmetry for which the classical motion 
is integrable and hence the Schrddinger equation is sep- 
arable. Symmetry-breaking terms are treated by pertur- 
bation theory. This approach fails when the symmetry- 
breaking terms become too large and typically numeri- 
cal solution is employed. An alternative approach which 
has been used successfully in atomic physics during the 
past decade is to use periodic orbit theory for the den- 
sity of states (DOS) of non-integrable systems. Until 
recently this approach has had few applications in con- 
densed matter physics and its power is not widely appre- 
ciated. Recently however Fromhold et al. |l| studied the 
tunneling I-V characteristic through a wide quantum well 
in a strong magnetic magnetic field oriented at an angle 
9 to the tunnel barriers. As they varied the voltage at 
fixed magnetic field (B = 11T) they found that resonant 
tunneling peaks (corresponding to the sub-band thresh- 
olds at 9 — 0) would abruptly change in frequency. Sim- 
ple estimates show that this behavior occurs well outside 
the regime in which sub-band mixing can be treated by 
perturbation theory. Modelling of the classical dynam- 
ics in this case indicated that the motion was completely 
chaotic and the new structure in the tunneling density 
of states (DOS) was explained by the appearance of new 
unstable periodic orbits via Gutzwiller's Trace formula 
0. A very important later experiment looked at the 
I-V peaks in the entire (plane) parameter space of mag- 
netic field and voltage for relatively small tilt angles || . 
At small tilt angles the frequency-doubling behavior was 
found at much lower magnetic field (B « 5T) but was 
"re-entrant" , i.e. would disappear at yet higher fields 
(Fig.l). At higher tilt angles a more complex pattern of 
peak-doubling and tripling is observed ||. This exper- 
iment was of particular interest because it mapped the 
transition regime to chaos in which stable and unstable 
periodic orbits exist and evolve according to the general 
constraints governing the transition to chaos. No previ- 
ous experiments of this type exist in condensed matter 



physics. 

Stimulated by this experiment Shepelyansky and Stone 
Q modelled its non-linear dynamics and showed that the 
system underwent a KAM transition to chaos described 
by a 2D effective map similar to the standard map || 
in which the degree of chaos was controlled by a di- 
mensionless parameter (3 = 2vqB/E ~ B/y/V (where 
£o = to*i?q/2 rj eV is the total injection energy). Hence 
at fixed tilt angle the degree of chaos is constant along a 
parabola V ~ B 2 //3 2 , but increases with increasing B at 
fixed V 0. The initial frequency-doubling island occurs 
well before the dynamics becomes fully chaotic and ref. 
noted that its onset appeared to be associated with a 
period-doubling bifurcation of a stable period-one orbit 
which surprisingly restabilizes for even higher (3 (chaos 
parameter) . In this work we develop for the first time an 
analytic theory of the periodic orbits controlling the main 
DOS oscillations in the transition region to chaos. The 
orbits responsible for the peak-doubling at small tilt an- 
gles are different from those previously considered Q and 
account for the re-entrant magnetotunneling behavior at 
small both qualitatively and quantitatively. At higher 
tilt angles these orbits remain important but others pre- 
viously considered 0] make significant contributions as 
well. 

The non-linear conductance of the well is proportional 
to the local DOS of the well in the vicinity of the left 
(emitter) barrier. Electrons tunnelling into the well at 
high voltages gain kinetic energy as they accelerate in the 
field and collide with the collector barrier. Over several 
collisions in the well the electron loses this energy by optic 
phonon emission. Therefore the tunneling resonances are 
substantially broadened and only are sensitive to struc- 
ture in the DOS on energy scales > h/T opt ~ hmeV. 
The Gutzwiller Trace formula relates the contribution to 
the DOS which is oscillatory with energy to a sum of 
contributions from each periodic orbit (PO) Q; the op- 
tic phonon broadening implies that only the few short- 
est PO's will give resolvable structure in the experiment. 
The Trace formula also implies that stable PO's typically 
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give larger contributions than do unstable, so that in a 
mainly chaotic phase space the biggest DOS oscillation 
comes from the few remaining stable PO's. 

The origin of chaos in the tilted well is as follows. At 
zero tilt angle (9 = 0) the acceleration along the elec- 
tric field E = Ez and the transverse cyclotron motion 
decouple and are integrable. Collisions with the barri- 
ers (y z — > — v z ) do not transfer energy between the cy- 
clotron and longitudinal motion. Once the B field is tilted 
(B = B cos 9z + sin Qy) the electron executes cyclotron 
motion around the B direction, with a superimposed drift 
velocity Vd = (E / B) s'mBx, and accelerates along B due 
to the component E • B = Ecos(9). This motion is still 
integrable. However now collisions with the barriers in 
general do mix the cyclotron and longitudinal energies 
e c ,£l and make the total dynamics non-integrable. The 
amount of energy exchange between those two degrees 
of freedom Ae depends sensitively on the phase of the 
cyclotron rotation at impact; for example if the velocity 
falls precisely in the x — z plane Ae = 0. Since the rate of 
change of the phase is u> c = eB/m*, we find (J] increas- 
ing chaos with increasing B. Thus the restabilization of 
period-one orbits at higher B requires explanation. 

With the choice A = (-Bycos(6) + B*sin(0),O,O), 
the system can be described by an effective Hamiltonian 
with two degrees of freedom 
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where U(z < 0) = 0, U(z > 0) = oo. If e > eV, 
the electron can retain enough longitudinal energy on 
collision at z = to collide again with the emitter bar- 
rier. But if £j = Eo — eV <C eV, (as is the case for the 
experiments of refs. [|]J^]), then |Ae| is much larger at 
the collector barrier (z = 0) than at the emitter barrier. 
Therefore the most salient features of the classical dy- 
namics can be described by the'T-barrier model" (IBM) 
[^j in which the emitter barrier U(—z — d) is omitted in 
Eq. (|]) above M. The periodic orbit structure of the 
IBM depends only on (3,9 [§. This can be seen directly 
from Eq. ([!]): if time is measured in units of uj^ 1 and 
distance in units E/Buj c the rescaled IBM Hamiltonian 
depends only on (3, 9. 

For 9 = the system is integrable and all the periodic 
orbits are of two types. First, degenerate families corre- 
sponding to lu c /ujl — p/q where p, q are integers. Here 
ujl = iteE I \J2mEh is the frequency of the periodic mo- 
tion in the longitudinal direction E = z, of the uniformly 
accelerated electron bouncing off the collector barrier. 
Note that this frequency decreases with increasing lon- 
gitudinal energy as 1/ ' y/eZ to its minimum value lul = 
2ttuj c /(3 when all the energy is longitudinal (el = eo)- Be- 
low this energy (determined by (3 m in — p/q) these helical 
orbits (HO) cease to exist. However there is one isolated 



periodic orbit which exists for all energies. This travers- 
ing orbit (TO) has no transverse energy (el = £o)j its 
period is To = 2mvo/eE = (3/uj c . It is is distinct from all 
helical orbits except if lo c Tq = f3 = 2np (p integer) when 
an infinitesimal increment of transverse energy generates 
a one-bounce HO of period To. The TO determines the 
sub-band energy spacings of the triangular well by (in- 
tegrable) semiclassical quantization, SE = Ti/Tq and will 
be relevant for the DOS oscillations at 9 ^ 0. 

First, all one-bounce orbits for 9 ^ must be "non- 
mixing" , i.e. there can be no exchange of El and e c at 
the bounce as these are separately conserved between 
bounces and the orbit is periodic. It follows that all 
such orbits must collide with with v y = 0, since this 
is the condition for zero energy exchange. This geo- 
metric property of all the one-bounce orbits allows us 
to to derive their periods, T, analytically. Consider 
such an orbit in the frame of reference moving with the 
drift velocity Vd = Esm(9)/B. In this frame the or- 
bit projected on the plane perpendicular to B is a por- 
tion of a circle with a chord removed of angular size 
20o = 2 arcsin(u ( iCi; c T/2w c ), where (j)o_ is the cyclotron 
phase just after the collision, and v c is the cyclotron 
velocity. Therefore, v c = {v d Lu c T/2) sin((27r - lu c T)/2). 
Since the motion along B between collisions is uniform 
acceleration, at the barrier |v • B| = eE cos(9)T/2m. Us- 
ing energy conservation we obtain (for arbitrary 9) 



(/3/2) 2 - (a; c T/2) 2 
(l-(w c T/2)cot(w c T/2)) : 



= sin 2 (9). 
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The solutions of Eq. (g|) are shown in Fig. 2. Note 
that for (3 < 2tt there exists only one orbit with pe- 
riod uj c To fa f3 as for the TO of the untilted system, 
and two more orbits appear at thresholds given by /?* ~ 
2itp, p = 1,2, . . .oo. Thus we have an infinite set of 
tangent bifurcations ]i| in which two new orbits (one sta- 
ble and one unstable) appear; for small 9 we find from 
Eq. (§) (31 = 2np{l + (9/irp) 1 / 3 + 0((9 4 / 3 )). It follows 
from Eq. (^) that for arbitrary (3 there exist 2N (with 
N = integer(/3//3^r)) such orbits whose periods are ap- 
proximately an integer multiple of 2tt/lo Ci and one orbit of 
longest period, To, which approaches the period co c Tq = (3 
of the unperturbed TO. This orbit is a recognizable de- 
formation of the TO at 9 = 0(see inset, Fig. 2a), whereas 
the other orbits are deformations of the 9 = helical or- 
bits of period oj c T = 2irp. The asymptotes at lo c T = 2irp 
arise because for 9^0 the period T = 2irp/oj c is not a 
solution of Eq. (||) for any /3. If it were a solution then 
in the drift frame the orbit would be a full circle, which 
is impossible due to the collision. The helical orbits with 
lu c Tq <C (3 have substantial cyclotron energy and typically 
remain far from the emitter barrier (see inset, Fig. 2a); 
hence the TO is the only period-one orbit important for 
the tunneling DOS and it defines the "sub-band" spacing 
of the tilted well fl. 
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The stability of these periodic orbits is described by 
a monodromy matrix || M; the orbit is unstable if 
|Tr(M)| > 2. For 1-bounce orbits we find 

Tr(Mi) = 4cos 4 (6>) (tan 2 (6») + (u c T/2) cot(w c T/2)) 

x (tan 2 (9) + sin(uj c T)/(uj c T)) - 2 (3) 

The results (plotted in Fig. 2) show that for small 
tilt angles and (3 < 4tt (the experimental parameter 
range) the TO goes unstable only in small intervals of 
(3 around (3 = tt, Ztt, .... This can be understood by con- 
tinuity. One can easily show that at 9 = the TO is 
only marginally stable (|Tr(M)| = 2) at /3 = lu c T = 
irp, p = 1,2... and is otherwise stable. Since the tran- 
sition to chaos is smooth jj] with increasing 9, instabil- 
ity can only appear for 9 <C 1 in small intervals around 
f3 = up. At (3 = f3p ~ 7r(2p-l)(l + 36» 2 /27r 2 (p-l/2) 2 the 
instability occurs by period-doubling bifurcation and a 
rapid restabilization occurs at (3? 1 ps n(2p — 1)(1 + 6> 2 (2 + 
1/-K 2 {p — 1/2) 2 ). At (3 ~ 2irp as noted above a stable 
and unstable period-one orbit appear so that two stable 
period-one orbits coexist; the one with the longer period 
is the TO which hence "jumps" in period. (The discon- 
tinuity in the peak data near B ss 10T, V ss 0.4V in Fig. 
1 occures right at [3 = (3\.) After this jump the helical 
orbits rapidly destabilize with increasing (3. For large 9 
or (3 the TO never restabilizes after its bifurcation near 
(3 ~ 7r, 37r, . . ., but it always reappears as a stable orbit at 
the tangent bifurcation near [3 ss 2ixp. Thus for arbitrary 
9 the TO undergoes an infinite sequence of destabiliza- 
tions and restabilizations and stable period-one islands 
exist at arbitrarily high chaos parameter. 

How does this explain the re-entrant behavior of the 
experimental DOS? As the magnetic field is increased (3 
first passes through (3 = (3 d « it and the TO bifurcates: 
for some interval of B there exists a stable period-two or- 
bit (T2 « 2To). The Trace formula implies that if a stable 
period-two orbit appears there should be a large increase 
in the amplitude of the DOS oscillation at Ti/2To (see 
Fig. 3) giving the period-halving. Eventually the period- 
two orbit destabilizes by period-doubling bifurcation, its 
contribution to the DOS is lost, and there is a return to 
the original sub-band spacing determined by the TO (see 
Fig. 1). However in contrast to simple area-preserving 
2D maps of the Chirikov type |^| the increase in chaos 
is not monotonic; a new stable TO appears by tangent 
bifurcation at yet higher B and its period-doubling bifur- 
cation defines the next interval of peak-doubling as seen 
in Fig.l. 

Although the IBM explains many of the puzzling qual- 
itative features of the data of ref. || , for further insight 
and quantitative comparison we now consider the two- 
barrier model, Eq. (0). The PO's in the 2BM depend 
not only on (3,9, but also on 7 = £o/eV. We calculate 
7 w 1.1 — 1.2 for the experiment of ref. [1 based on 
the standard model for the emitter state E ] . Using the 



non-mixing property of the TO we can again calculate its 
period and stability analytically jll| and (assuming con- 
stant 7) we plot the result in Fig. 2b vs. (3. Note that in 
the 2BM for 9 = the TO obeys uj c T = (3(1-^1 - 1/7) 
(not (3) and the behavior for small 9 approaches this line 
except near (3 = 2np where the helical period-one or- 
bits appear. These orbits have too little energy to reach 
the emitter barrier and hence are exactly as in the IBM. 
The TO also has a very similar geometry as in the IBM 
except of course for most (3 it collides with both barri- 
ers. However a significant difference arises in its stability 
properties. The order of the bifurcation and restabiliza- 
tion points is interchanged as a function of (3 (or B) (see 
insets, Fig. 2a,2b). This implies that the period-doubling 
region occurs for (3 < (3 d (smaller B) and not for (3 > (3 d 
as in the IBM. This difference can be understood by not- 
ing that the period-two orbit requires substantial e c and 
hence occurs only above a threshold energy, £2, deter- 
mined by Tq{e2,B) ps tt/lj c . In contrast to the IBM, 
in the 2BM To is a decreasing function of energy since 
a higher energy collision with the emitter barrier just 
returns the electron to the collector sooner. Thus the 
high-energy region in which a stable period-two orbit ex- 
ists, occurs when lu c Tq < tt, i.e. for smaller values of [3. 
Denoting the point of bifurcation by /?2 , the restabiliza- 
tion occurs for [3r < @2- At /3 3 < (3r the period- two orbit 
destabilizes and at (3± < P3 it disappears in an (inverse) 
tangent bifurcation with a period-two orbit which doesn't 
reach the emitter. Thus we identify (3^ with the onset of 
period- halving in the tunneling oscillations in Fig. 1, and 
/?2 with the return to period-one oscillations. The same 
sequence occurs again at higher [3 = f3' 2 , f3' 4 . These peak- 
doubling boundaries are calculated analytically from the 
2BM using the estimated experimental parameters (with 
no fitting) and are found to agree rather well with exper- 
iment (Fig.l). 

The peak behavior at higher tilt angles is more com- 
plicated, as the TO is unstable over a large interval of 
[3 and the contribution of it and its offspring to the tun- 
neling DOS needs to be compared with other unstable 
orbits. Comparison to the data of ref. Q at 38° shows 
peak-doubling and tripling regions in very good agree- 
ment with the bifurcation and trifurcation of the TO 
and a re-entrant single-peak region at high B in reason- 
able agreement with the restabilization of the TO after 
tangent bifurcation |ll[] . Other regions of multiple peaks 
appear to be associated with previously identified orbits 
JlJ which are not bifurcations or trifurcations of the TO 
but whose interval of existence can be understood in a 
similar framework [ Q |. 

In Fig. 3 we confirm that the DOS for the tilted well 
does show peak-doubling in the interval predicted by our 
theory. Here we show a numerical calculation of the oscil- 
latory part of the DOS appropriately scaled using Eq. (1) 
pr| for 9 = 20° . As the energy increases through the clas- 
sical bifurcation point the DOS oscillations clearly pick 
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up a secondary maximum corresponding to the period- 
halving, but at even higher energies the secondary max- 
ima disappear at the predicted restabilization point for 
the TO. Note in Fig. 1 that the secondary maxima are 
not observed at lower and higher B, V where the classical 
theory implies a stable period-two orbit exists. However 
at lower fields the period T2 ~ l/^ c increases, so the op- 
tic phonon damping ~ exp(— T2/T op t) removes the contri- 
bution of the period-2 orbit from the tunneling spectra. 
Hence secondary maxima die out in the period-two layer 
at low fields (see lower curve of Fig. 3). At higher B. V 
we propose a different mechanism suppresses the peak- 
doubling. Typically, the period-two orbit has a substan- 
tial amount of cyclotron energy e c ~ eo — eV = (7— l)eV, 
which is greater then the cyclotron energy in the emitter 
state (~ Ticu c ). Tunneling to the period-2 orbit therefore 
requires energy transfer from longitudinal to cyclotron 
motion. The amount of this energy transfer increases 
with increase of B and V, leading to a decrease in the 
tunneling matrix element for this process at high fields 
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FIG. 1. Resonant tunneling peak locations for 9 = 11° 
from [3]. Pairs of solid curves enclose the regions of exis- 
tence of period-two orbit for 2BM. Dashed line corresponds 
to the destabilization of the period-2 orbit. Variation of the 
effective mass with energy has been taken into account. 
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FIG. 2. (a) Periods of one-bounce orbits vs. /3 = 2{vqB/E) 
for one-barrier model with 9 = 11°. Solid lines denote stable 
orbit, dotted unstable. Dashed line is 9 = result for TO. 
Triangles = period-doubling bifurcations, circles = restabi- 
lizations. Insets: (top left): Traversing orbit (TO) for j3 = 10 
projected in y-z (left) and x-y (right) plane, (top right) same 
for helical orbits (HO). (Lower right): blow-up near j3 — it. 
(b) Same for two-barrier model with e/eV — 7 = 1.1. Lower 
dashed line is 9 = result for 2BM, upper for IBM. Grey-scale 
regions denote intervals of existence of period-two orbit which 
is born in tangent bifurcation and disappears in a backwards 
period-doubling bifurcation. 
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FIG. 3. Oscillatory contribution to DOS calculated numer- 
ically vs. action So(eo)/h of unperturbed TO. Note onset of 
secondary maxima at classical bifurcation energy (1st arrow) 
and disappearance at classical restabilization energy (2nd ar- 
row). (Bottom Trace, B = 3.53T, V = 0.2V, 6 = 20° ): stable 
period-two interval has decreased and shifted to lower ener- 
gies and secondary maxima are no longer apparent (see text). 
Dashed line is ratio of the period-two to period-one contri- 
butions to the DOS for B = 5T, V = 0.4V, 6 = 20° from 
Eqs. (2), (3) using the Trace formula [2]. Level-broadening is 
h/Topt = 6.6meV. 
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